Renormalization group analysis of the M-p-spin glass model with p = 3 and M = 3 
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We study an Af-p-spin spin glass model with p — 3 and M = 3 in three dimensions using the 
Migdal-Kadanoff renormalization group approximation (MKA) . In this version of the p-spin model, 
there are three (Af = 3) Ising spins on each site. At mean-field level, this model is known to have 
two transitions; a dynamical transition and a thermodynamic one at a lower temperature. The 
dynamical transition is similar to the mode-coupling transition in glasses, while the thermodynamic 
transition possibly describes what happens at the Kauzmann temperature. We find that all the 
coupling constants in the model flow under the MKA to the high-temperature sink implying that 
the mean-field features disappear in three dimensions and that there is no transition in this model. 
The behavior of the coupling constant flow is qualitatively similar to that of the model with p = 3 
and M — 2, for which only a single transition is predicted at the mean-field level. We conclude 
that for p-spin models in three dimensions, fluctuation effects completely remove all traces of their 
mean-field behavior. 



Understanding the dynamics and thermodynamics of 
supercooled liquids and structural glasses is one of the 
most chahenging problems in condensed matter physics. 
Theoretical advances have been made using an idea ob- 
tained from certain mean- field spin glass models [ij . One 
example is the infinite-range p-spin spin glass model [l| , 
which shows two phase transitions. There is a dynam- 
ical transition at a temperature below which an er- 
godicity breaking occurs. The connection to structural 
glasses comes from the similarity of the dynamical equa- 
tions of these spin glass models near Td to those in the 
mode-coupling theory of the structural glass transition 
0. Below the dynamical transition, a thermodynamic 
transition takes place at Tq with the low temperature 
phase described by the one-step replica symmetry break- 
ing (RSB). This transition has been associated with a 
possible transition at the Kauzmann temperature Tr- Q 
in structural glasses. An important question is then how 
this mean-field picture gets changed in finite dimensions. 
Beyond the mean-field level, the dynamical transition at 
Td is expected to disappear because of activation pro- 
cesses over finite energy barriers separating metastable 
states. However, the thermodynamic glass transition Tq 
could still exist in principle in finite dimensional systems. 
In this paper, we use the Migdal-KadanofF (MK) real 
space renormalization group (RG) method to investigate 
this issue. 

The original p-spin glass model mentioned above, how- 
ever, is inconvenient to study on a simple hypercubic 
lattice. In Ref. Q, a new class of p-spin glass mod- 
els called the M-p-spin glass model has been proposed. 
These models consist of Ising spins on a hypercubic lat- 
tice with nearest-neighbor p-body interactions. The new 
feature of these models is that there are M Ising spins 
on each site. As the value of the parameters M and p 
change, these models exhibit two distinct types of behav- 
ior at mean-field (MF) level. For example, when p = 3 



and M = 2, the model at mean- field level has a single 
continuous transition to a full RSB state. On the other 
hand, for p = 3 and Af — 3, one has the two transi- 
tions described above. One can actually classify all these 
models for general p and M , into the two mean- field tran- 
sition types [1]. Therefore the M-p-spin glass model is 
an ideal testbed for studying how thermal fluctuations 
change mean-field results. The p = 3, M = 2 model has 
been studied by Monte Carlo simulations in four dimen- 
sions 0] and the results were interpreted as evidence for 
a single continuous transition. The Migdal-Kadanoff RG 
has been applied [^, 0] to this model, using versions to 
mimic both three and four dimensions, and it was shown 
that the coupling constants flow to the high-temperature 
sink implying that there is no transition in both these 
dimensions. This model has been shown Q to be in the 
same universality class as the Ising spin glass in a field, 
which is consistent with the absence of a phase transition 
in less than six dimensions Q . The correlation length ob- 
tained from the Migdal-Kadanoff RG grows quite rapidly 
with decreasing temperature to a value of order 30 lat- 
tice spacings. Thus in a numerical simulation of a finite 
system whose linear extent is of order 10 lattice spacings, 
it is therefore very likely that the rapidly growing corre- 
lation length seen would appear to be implying there is 
a transition but this would just be an artifact of fi- 
nite size effects. The non-mean-field regime of the p — 3 
and M — 2 model has also been simulated using a one- 
dimensional long-ranged Ising spin glass where the cou- 
plings decrease with their separation rij as ^/^ij- By 
varying the exponent a one can effectively tune the one- 
dimensional system to have similarities to the short-range 
problem in spatial dimension d [9]. The mean- field tran- 
sition was shown to be destroyed at values of a which 
correspond to low values of d. 

It is a purpose of this paper to find whether the mean- 
field transitions in the p = 3 and M = 3 model, which 



are relevant to discussions on structural glasses, survive 
to three dimensions using the Migdal-KadanofF real space 
RG. The MKA is an approximate RG procedure which 
works best in low-dimensions d such as two and three 
where it gives excellent results for the Ising (i.e. p = 2) 
spin glass. 

The M-p-spin model we study in this paper is charac- 
terized by M Ising spins (Tj-"\ a = 1, 2, • • • , M on each 
site I of a hypercubic lattice. The Hamiltonian is given 
in terms of products of p spins chosen from the spins in a 
pair of nearest neighbor sites. In this paper, we focus on 
the p = 3 case and the Hamiltonian of the 3-spin glass 
model is given by 

M M 
(ij) a<l3 7 

(1) 

where the notation {ij) means that the sum is over all 
nearest neighbor pairs i and j. Note that the number of 
different coupling constants, J-"^^'^ and J^^"^'' for given 

{ij) is just 2M(*^) = M2(M - 1). All these couplings 
are chosen independently from a Gaussian distribution 
with zero mean and width J. Therefore for the 3-spin 
glass models, there are 4 and 18 independent couplings 
for M — 2 and M = 3, respectively. 

The above Af-p-spin model can be put into a field the- 
oretical framework. A standard way to do this is to use 
the Hubbard- Stranotovich transformation on the repli- 
cated partition function and then trace over the spins. 
The resulting field theory associated with this model is 
the following Ginzburg-Landau- Wilson Hamiltonian 

^GLW - /rf'r|l^[Vg„,(rr + i^g^,(r) 

-^Tr,3(r)-!|5],3^(r)|, (2) 

a<b ) 

where Qofclr) is the order parameter and a and b are 
replica indices running from 1 to n with n — > 0. At mean 
field level, this model is known to show very different 
behavior depending on the value of i? = W2/W1. When 
R > I, there are two transitions at the mean field level 
as described above; a dynamical transition at some tem- 
perature Td and a thermodynamic transition at a lower 
temperature To to a state with one-step replica symme- 
try breaking. The two transitions have been discussed 
extensively in connection with what happens in struc- 
tural glasses near the mode-coupling and the Kauzmann 
temperatures, respectively. For the case where i? < 1, 
only a single transition to a state with full RSB is ex- 
pected at mean- field level. In Ref. [1], the ratio R was 
evaluated for the M-p-spin model for general values of 
M and p. The cases we are interested in this paper. 




n 



FIG. 1: RG flow of the three-spin coupling strength 
K{n)/K{0) normalized with its initial value at iteration step 
n. The upper panel is for p = 3 and M — 3, and the lower one 
for p = 3, Af = 2. One each panel, the curves correspond from 
bottom to top to temperatures T/ J — 10, 8, 6, 4, 2, 1, 0.5, 0.25 
and 0.1. 

namely p = 3, A/ = 2 and p = 3, M = 3, correspond 
to i? = 3(1 - 1/V2) ~ 0.879 and i? = 2, respectively. 
Therefore the two models show very different mean-field 
behavior. It is the purpose of this paper to investigate 
how thermal fluctuations affect the two kinds of mean- 
field behavior in three dimensions. 

To do this we apply the Migdal-Kadanoff real space 
RG to these models. We follow closely the approximate 
bond moving scheme described in detail in Ref. For 
a three dimensional cubic lattice, four bonds are put to- 
gether to form a new bond. The coupling constant of this 
new bond is just given by the sum of the coupling con- 
stants of the four bonds. We can obtain a coarse grained 
lattice by taking a trace over the spins at the site con- 
necting two new bonds. The decimation procedure can 
be continued n times for a system of size L = 2". Al- 
though we initially start from only the three-spin cou- 
plings K^"/^)-^ = j("Ph/T and K'^(°'P^ = Jii^'P) /T 
drawn from a Gaussian distribution with standard devia- 
tion (width) K = J /T at temperature T, the decimation 
procedure generates additional types of couplings as well 
as on-site "field" terms which involve the spins only on 
one site. In the Migdal-Kadanoff analysis, we have to 
keep track of all these couplings. 

For the p = 3, Af = 2 case, on each bond there are ini- 
tially four three-spin couplings. In addition to those, we 
have to consider four two-spin couplings and one coupling 
that involves all the four spins on the nearest-neighbor 
sites. There are also four field terms on the two nearest 
neighbor sites and two on-site terms involving two spins 
on one site. If we include a constant contribution to the 
Hamiltonian, the renormalized Hamiltonian carries a to- 
tal of 16 coupling constants. For the p = 3, M = "i case, 
on the other hand, each bond is connected to 6 spins, 3 
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FIG. 2: RG flow of the two-spin coupling C(n) and the four- 
spin coupling D{n) of the p = 3, Af = 3 model at the n-th 
iteration step. The curves correspond from bottom to top to 
temperatures T/J = 8, 4, 2, 1, 0.5 and 0.25. 



spins on each site. There are couphngs involving 2, 3, 4, 
5 and 6 spins from a pair of nearest neighbor sites, whose 
standard deviations are denoted by C, K, D, E and F, 
respectively. The numbers of independent coupling con- 
stants connecting 2, 3, 4, 5 and 6 spins are given by 9, 18, 
15, 6 and 1, respectively. There are also on-site terms. 
The numbers of on-site fields involving 1, 2 and 3 spins 
are 6, 6 and 2, respectively. Including a constant term, 
we have to keep track of a total 64 couplings in the case 
oi p — i and M — 3. 

A key step in the calculation is setting up and solv- 
ing the recursion relation resulting from the decimation 
procedure. For p = 3 and M = 3, this amounts to solv- 
ing a system of 64 simultaneous equations for a new set 
of couplings, which requires an inversion of a constant 
64 X 64 matrix. This only needs to be done once before- 
hand and not for every iteration. In the actual numerical 
calculations, 10^ bonds are prepared. On each bond 18 
three-spin couplings, K^°'^'>'' and K'^^°'^'> are chosen inde- 
pendently from the Gaussian distribution with zero mean 
and width K{0) ~ J /T aX temperature T. All the other 
couplings including the fields are set to zero initially. We 
then randomly choose two sets of 4 bonds to form two 
new bonds. Using the recursion relation, we obtain a 
set of 64 renormalized couplings. This procedure is con- 
tinued until we get 10^ new bonds, which completes the 
first iteration. As we iterate the same procedure, we can 
obtain the flow of 64 couplings. At each step of the itera- 
tion, we evaluate the width of each distribution. We note 
that there is a certain freedom in where to move on-site 
flelds when the four bonds are put together. The main 
results do not depend on the method chosen. We follow 
the prescription given in Ref. when three bonds are 
moved to combine with a bond, the fields on the three 
bonds are placed at the site that is to be traced over. 



The RG flow of the width of the three-spin coupling 
K{n) at the n-th iteration normalized with the initial 
value K{Q) = J/T is shown in Fig. [1] for the 3-spin glass 
model with M = 3 and M = 2, respectively. Somewhat 
surprisingly, both flgures show essentially the same fea- 
tures. At all temperatures, the couplings decrease to zero 
in the long length scale limit indicating an RG flow to the 
high temperature sink. This implies that there is no flnite 
temperature phase transition in these models in contrast 
to the corresponding mean-field pictures. At high tem- 
peratures, the coupling simply decreases to zero. As the 
temperature is reduced, however, the coupling increases 
during the first few iterations reaching a maximum at 
the second or third iteration, but eventually decaying to 
zero. At very low temperatures T/J < 0.25, the flow 
patterns of K{n)/K{0) reach an asymptotic form. The 
only difference between M — 3 and M — 2 cases is that 
the couplings in the AI — 3 model decay to zero more 
slowly as the iteration proceeds. We can deduce from 
this that the correlation length in the M — 3 model is 
longer than that of the M = 2 model. 

The other couplings for the p — 3 and M — 3 model 
also flow to zero as can be seen from Fig. [2] The upper 
panel shows the RG flow of the width C(n) of the two- 
spin couplings for various temperatures, which reaches 
a maximum at the first iteration and then decreases to 
zero as the iteration proceeds. The four-spin coupling 
strength D{n) also shows the same behavior except that 
it has the largest values at the second iteration step. The 
five- and six-spin couplings (not shown here) also exhibit 
essentially the same behavior. The on-site fields, how- 
ever, do not decrease but increase and saturate to a con- 
stant value in the long length scale limit as shown in 
Fig. m The same behavior was observed in the Migdal- 
Kadanoff analysis of the p — 3, M = 2 model. Therefore, 
the p-spin glass model for both M — 2 and AI — 3 cases 
correspond on large length scales to a model with random 
fields acting on spins on decoupled sites. 

We now study how the standard deviations of cou- 
plings behave as a function of the length scale L = 2". 
The standard deviation K of the three-spin couplings 
shown in Fig. [T] decays as exp(— L/^(T)), where ^(T) 
at temperature T can be interpreted as the correlation 
length. In Fig. IH the correlation lengths of the two 3- 
spin glass models with M = 2 and AI — 3 are shown. 
The correlation length for the p = 3, M — 2 model has 
been already obtained in Ref. 0]. The AI = 3 model 
has longer correlation length than the AI — 2 model at 
any given temperature. The correlation lengths of the 
two models basically show the same temperature depen- 
dence. They increase quite rapidly as the temperature 
is lowered and then saturate at large values in the zero 
temperature limit. 

In summary, we have investigated the RG flows of 
the coupling constants in the M-p-spin glass model with 
p = 3 and M = 3 in three dimensions using the Migdal- 
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FIG. 3: RG flow of on-site fields at temperature T / J = 0.5 
for p = 3 and M = 3. Circles, squares and diamonds corre- 
spond to the one-spin, two-spin and three-spin one-site terms, 
res 




FIG. 4: Correlation length ^ as a function of temperature 
T / J for the p = 3, M = 3 model (filled circles) and the p = 
3, M = 2 model (open circles). 

KadanofF approximation. All the coupling constants flow 
to the high temperature sink implying that there is no fi- 
nite temperature phase transition in this model. The cor- 
relation length increases quite rapidly as the temperature 
is lowered and saturates in the zero temperature limit. 
All these features are very similar to those of its AI = 2 
counterpart. The only difference is that the M ~ 'd model 
has longer correlation lengths. This is quite remarkable 
in that the two models show completely different behav- 
ior at the mean-field level. Therefore we conclude that in 
the 3-spin spin glass models, thermal fluctuations in three 
dimensions completely destroy the transitions predicted 
by mean-field theory: whether i? > 1 or is less than 1 
results in only a quantitative change in the magnitude of 



the correlation length. 

What stimulated this investigation was the renormal- 
ization group analysis by Cammarota et al. of a 

replicated Hamiltonian similar to Hglw in Eq. Like 
us they used a MK approximation appropriate to three 
dimensions, but instead found a transition which they 
identified as the transition expected in random first-order 
transition theory pj at a temperature they identified as 
Tk- We believe the difference in results stems from their 
decimation step in the MK procedure. In our work it 
is carried out exactly. In the calculation of Ref. 
however, it was done by a steepest-descent approxima- 
tion to an integral over n x n replica matrices, and the 
magnitude of the corrections to this approximation were 
not assessed. They themselves suggested that the cal- 
culations actually done in this paper should be carried 
out as a check of their decimation approximation. Our 
calculations alas provide no support for it. 
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